Charge symmetry breaking in np — > dir° 
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The forward-backward asymmetry in np — + d% , which must be zero in the center-of-mass system 
if charge symmetry is respected, has been measured to be [17.2 ± 8(stat) ± 5.5(sys)] x 10~ 4 , at an 
incident neutron energy of 279.5 MeV. This charge symmetry breaking observable was extracted 
by fitting the data with GEANT-based simulations and is compared to recent chiral effective field 
theory calculations, with implications regarding the value of the u d quark mass difference. 



PACS numbers: 11.30.Er, 13.75.Cs, 24.80.+y 

In the quark model, the breaking of charge indepen- 
dence and charge symmetry arises from the mass differ- 
ence of the up and down current quarks and the elec- 
tromagnetic interaction between quarks. The basic np 
interaction is particularly sensitive to such fundamental 
effects since the "background" Coulomb force is absent 
in this system. Indeed, charge symmetry breaking (CSB) 
has been unambiguously observed 0, 0- 13 m n P elas- 
tic scattering at three different energies. Measurement 
of CSB in the inelastic np — > dir° reaction complements 
the existing data in that it is sensitive to contributions 
that are absent in the elastic channel. Furthermore, 
this reaction is unique as a testing ground for effective 
field theory calculations addressing the important issue 
of isospin symmetry violation in pion-nucleon scattering. 
The observable of interest in np — ► dir° is the center-of- 
mass forward-backward asymmetry, Ar, , which we define 
as 



A^e)^ 6) -^- 



tr(0) + a(ic - 6) ^ 

where 9 is the angle between the incident beam and the 
scattered deuteron. Note that the asymmetry must be 
zero if charge symmetry is conserved. We report on a 
measurement of this asymmetry at a neutron energy a 
few MeV above the reaction threshold (275.06 MeV) .and 
compare our result to recent theoretical predictions [J, |5j 
bearing on such fundamental questions as the u d quark 
mass difference and our understanding of QCD dynamics 
and symmetries in low-energy hadronic interactions. 

THE EXPERIMENT 

The experiment was performed at TRIUMF with a 
279.5 MeV neutron beam, a liquid hydrogen target, and 



the SASP magnetic spectrometer |6j positioned at 0°. 
With these near threshold kinematics and the large ac- 
ceptance of SASP, the full deuteron distribution from 
np — > dir° was detected in one setting of the spectrom- 
eter thereby eliminating many systematic uncertainties. 
These deuterons form a distinct kinematic locus in mo- 
mentum vs laboratory scattering angle, which is shown 
in fig. \I\for the collected data. 
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FIG. 1: Kinematic locus of np — » dir° data. 



The TRIUMF CHARGEX facility produced the 
neutron beam by passing a high intensity proton beam 
through a thin 7 Li target. A sweeping magnet deflected 
the primary proton beam into a well- shielded dump. The 
liquid hydrogen target (LH2) was centered 92 cm down- 
stream from the 7 Li target and was contained within a 



2 



flat cylindrical volume, 10 cm in diameter with a nomi- 
nal thickness of 2 cm. Two sets of veto counters (FEV1, 
FEV2) and a trigger counter set (FET) were each com- 
posed of a pair of plastic scintillators positioned above 
one another. This allowed more stable operation in the 
high (few MHz) particle rate environment. The thick 
veto scintillators were upstream of the LH2 and shad- 
owed it. The FET counters were positioned immediately 
downstream of the LH2. 

Three multi-wire proportional chambers, positioned 
upstream of the SASP entrance (FECs i.e. Front-End 
Chambers), provided tracking information for charged 
particles. Each FEC consisted of a pair of orthogonal 
wire planes. The first and last FECs, separated by 33 
cm, were mounted to measure vertical and horizontal co- 
ordinates. The third FEC was positioned midway be- 
tween the other two and rotated 40° with respect to 
them for efficiency measurements and to aid in multi-hit 
track reconstruction. Particle tracking near the SASP 
focal plane was provided by two vertical-drift chambers 
(VDCs). Three sets of scintillators, downstream from 
the VDCs, provided timing and particle ID information 
as well as sufficient redundancy to determine the efficien- 
cies of all focal plane area detectors. 

Measurements of np elastic scattering with incident 
neutron beams that filled the same target space and pro- 
duced protons that spanned the momentum distribution 
of the np — > dir° reaction provided a stringent test of the 
description of the spectrometer acceptance. Further de- 
tails on the apparatus and other technical aspects of the 
measurement are found in reference Ig. 



angle (as in fig. ^) and compared to a model which repre- 
sented the background due to C(n, d) reactions as a low- 
order polynomial and generated the locus of H(n, d)7r° 
events by Monte Carlo simulation of the beam, target, 
reaction cross section, spectrometer and detectors. 

The simulation was based on GEANT3. It began 
with a proton beam incident on the 7 Li target and in- 
cluded energy loss by the proton beam as well as the 
angular and energy distribution of neutrons from the 
7 Li(p, n) reaction. Production of deuterons according to 
the distribution of equation [21 was allowed in the LH 2 
target and other hydrogenous material such as scintil- 
lators and their wrapping. Standard GEANT tracking 
options were adopted for deuteron energy loss and mul- 
tiple scattering but the reaction losses, which amount to 
1-2% and are momentum dependent, were parametrized 
from data on deuteron elastic and reaction cross-sections 
from hydrogen and carbon 10]. Tracking through the 
SASP dipole used a field map obtained at 875 Amp and 
scaled up to the operating current of 905 Amp. Data 
were acquired in 10 different periods spanning two years 
and the simulation accounted for measured detector effi- 
ciencies, scintillator thresholds, missing FEC wires, and 
known changes in target thickness in a manner consistent 
with the actual running periods. 

To reduce the possibility of psychological bias in 
matching simulation to data, a blind analysis technique 
was used which incorporated a hidden offset to the A\ /Aq 
asymmetry parameter of the np — > dir° generator. The 
collaborators developing the simulation and extracting 
the observable did not know the value of the offset until 
all consistency checks had been satisfied. 



EXTRACTION OF An, 

Close to threshold, the np — > dir° cross-section in the 
center-of-mass frame is given by 

^-(0)=A o + A 1 P 1 (co S 6)+A 2 P 2 (cos6), (2) 
dil 

where Pi and P2 are Legendre polynomials. The A and 
A2 coefficients were previously measured at a number 
of energies within 10 MeV above threshold. The presence 
of charge symmetry breaking is reflected in the A\ term 
as it is odd in cos 9. In this standard parametrization, the 
angle integrated form of Afb is given by Aft, — kAi/Ao. 

For a given beam energy, cos 9 varies linearly with the 
longitudinal component of deuteron momentum in the 
laboratory reference frame. Ideally, the cos 9 distribu- 
tion would be found by a suitable, simple projection of 
the data of fig-H However, the measured deuteron locus 
is distorted by energy loss, multiple scattering, energy 
spread of the beam, and spectrometer acceptance mak- 
ing a direct extraction of Ar, impossible. Instead, the 
data were binned according to laboratory momentum and 



SYSTEMATIC EFFECTS 

The acceptance of SASP is a function of the initial 
target position and direction of the deuteron as well as 
its momentum. Non-uniformities in the momentum ac- 
ceptance of SASP would systematically produce a false 
asymmetry and had to be limited. High-statistics data 
from np elastic scattering were collected and compared 
to model simulations to determine a fiducial volume 
of uniform acceptance. For these calibration measure- 
ments, the SASP magnets were set to their values for 
the np — > dir° running, but the primary beam energy 
was adjusted so that the elastically scattered protons had 
a momentum deviation 6 = (p — Po)/po = —4, or +4% 
compared to the central momentum of the deuterons of 
interest. Projections of the np elastic data direction for 
position slices were formed, and the ratios of yields at 
—4% vs +4% were formed for both data and simulation; 
see fig. |5J The analysis software acceptance cuts in posi- 
tion and direction were then limited to the regions com- 
mon to both data and simulation which were uniform in 
momentum to the statistical precision of the data. 
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Simulation vs simulation comparisons were carried out 
to determine how strongly experimental parameters and 
other effects were correlated with Ai/Aq. For exam- 
ple, momentum dependent deuteron reaction losses and 
detection efficiencies are obvious mechanisms which can 
mimic the effect of a non-zero Ai/Aq. Combining each 
correlation with the independently-determined uncer- 
tainty of its parameter gave the systematic contributions 
shown in Tabled However, for the LH2 target thickness, 
the proton beam energy (Tb ca m) and the central momen- 
tum of SASP (po) the independent information was not a 
sufficient constraint. Therefore, these three parameters, 
along with Ai/Aq, were treated as free parameters and 
their values extracted from fitting the data. To this end, 
simulations were made and x 2 calculated for 81 points 
in a four-dimensional space, in which each of the four 
free parameters was stepped above and below a nominal 
value, x 2 minimization techniques 1 1 1 1 | were then used to 
obtain the values of the parameters at the global x 2 min- 
imum, while the local curvature of the x 2 surface gave 
their errors and mutual correlations. 
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FIG. 2: The vertical projection of the lab scattering angle for 
the center acceptance slice, elastic np scattering; A) data yield 
for S = -4%, no cuts; B) data ratio (— 4%/+4%), no cuts; C) 
data ratio (— 4%/+4%), full acceptance cuts; D) simulation 
ratio (— 4%/+4%), full acceptance cuts; E) ratio of C to D, 
normalized to the center bin. 



TABLE I: Systematic Error Contributions to An, 
Uncert x(l(T 4 ) 



FEV threshold 


2.5 


Separation between front and rear FECs 


2.5 


Longitudinal position of 7 Li 


2.5 


A 2 /A 


2 


Deuteron reaction losses 


1.5 


Detection efficiencies 


1.5 


Primary beam energy spread 


1 


Neutron angle 


1 


Background 


1 


FET threshold 


0.5 


Total 


5.5 



RESULTS AND DISCUSSION 

As a test of the model, the x 2 calculations and fitting 
were repeated on subsets of the data and simulated data, 
selected according to whether the reaction occurred in 
the Top or Bottom part of the LH2 target. A second 
test divided events into those originating in the Left or 
Right part of the target. The best fit values and errors of 
Ai/Aq (after removal of the offset) and the other three 
parameters are presented in Table ITT1 

The root mean square (rms) systematic error for the 
full acceptance and the four subspaces is ~ 2.7% with 
the standard binning scheme of 50 bins in 5 and 20 bins 
in 0, indicating a substantial discrepancy between data 
and the simulation. Pixel by pixel examination of the 
contribution to x 2 revealed a systematic difference in the 
profile of the locus along lines of steepest ascent. The sign 
of the differences tended to be positive at the peak and 
negative at both the "inner" and "outer" margins of the 
locus, possibly due to inadequate treatment of deuteron 
scattering in the simulation. 

A change in Ai/Aq will not change the ratio of counts 
in peak vs margins of the locus because it multiplies 
cos(# cm ). In contrast, the LH2 thickness, pq, and Tbcam 
all shift or broaden the locus and thus are sensitive to the 
ratio of locus counts at the peak vs margins. It is rea- 
sonable to expect further rebinning to remove sensitivity 
to unimportant details of the simulation without losing 
sensitivity to Ai/A . We repeated the x 2 grid search us- 
ing 20 bins in S and 10 bins in 8, and again with 10 bins 
in S and 5 bins in 9. As expected, the fractional error 
dropped to 2.1% and 1.4%, respectively, with Ai/Aq re- 
maining consistent within errors. A more sophisticated 
binning scheme which treated the locus as a set of "ellip- 
tical" and "radial" bins on top of rectangular background 
bins produced an rms error of 2.1% for 2500 background 
bins and 36 locus bins, and an rms error of 1.2% for 100 
background bins and 6 locus bins. In all fits and binning 
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TABLE II: Stability of the four free parameters over target 
subspaces; b = bottom; t = top; 1 = left; r = right. 



(Ai/Ao) relative relative 
(1(T 4 ) LH 2 (mm) p (MeV/c) 



relative 
T beam (MeV) 



full 34.4 ± 16 

b 30 ± 26 

t 20 ± 20 

1 29 ± 23 

r 15 ± 22 



0.94 ± 0.05 0.365 ± 0.015 0.048 ± 0.001 
0.39 ± 0.09 0.547 ± 0.025 0.086 ± 0.002 
1.14 ± 0.07 0.236 ± 0.018 0.021 ± 0.002 



1.21 ± 0.08 
0.75 ± 0.08 



0.273 ± 0.021 
0.427 ± 0.021 



0.042 ± 0.002 
0.051 ± 0.002 



schemes the best fit values of the asymmetry in the accep- 
tance subspaces agreed within errors with the value for 
the full acceptance, which is (34.4 ± 16) x 10~ 4 , implying 
Afb = [17.2 ± 8(stat) ± 5.5(sys)] x 1CT 4 . 

Theoretical predictions of Aft, have been made by 
Niskanen 4] using a meson-exchange coupled-channel 



model which showed that the major contribution by far 
is due to ttt] (and in]') mixing in both the exchange and 
produced (outgoing) meson. At our energy the prediction 
is Aft, = —28 x 10~ 4 , when accepted values are used for 
the rjNN coupling constant and the ir°r] mixing matrix 
element (<7^jvAr/4 7r = 3-68 from meson exchange NN po- 
tential models [13 and (n°]H\r]) = -0.0059 GeV 2 from 
analysis of r\ decay data |l3j). More recently, A^ was 
revisited []| within the framework of chiral effective field 
theory where the issue of charge symmetery breaking in 
the rescattering amplitude of the exchanged pion was ad- 
dressed. The resulting additional contribution to An, is 
then expressed in terms of two parameters 5m n and 5m^ 
representing contributions from the u d quark mass differ- 
ence and from electromagnetic effects within the nucleon, 
respectively. Specifically, at our energy, Afb is expressed 
as 



J 



ifb 



-0.28% x 



9r)NN 



t/47t(3.68) 



-0.0059 GeV 2 



0.87 

MeV 



(SrriM 



5m 



N 



(3) 



where the first term arises from 7rr/ mixing and the sec- 
ond from 7T° N scattering. With the introduction of the 
new term, Ar, changes sign and becomes positive with 
an estimated upper value around +69 x 10~ 4 , when large 
but reasonable values of 8m ^ and 5m^ are used @ . Our 
positive experimental result strongly suggests, therefore, 
that such isospin violating tt°N interactions as outlined 
in reference [5j are indeed significant. 

The parameters Sm^ and 5m n are also constrained by 
the proton-neutron mass difference as 



from the National Research Council of Canada. 



A w = 



m„ — m, 



= 5m N + 5m N = 1.29 MeV. 



(4) 



When our Ar, result is combined with equations [3] and 
0] and the values given above for the ijNN coupling 
constant and tttj mixing matrix element, we find that 
5m N = 1.66 ± 0.27 MeV and 5m N = -0.36 ± 0.27 MeV, 
assuming no theoretical uncertainties. We emphasize, 
however, that this last exercise is only meant to illustrate 
the significance and potential important implications of 
our Afb result. Further theoretical studies are currently 
underway |14l | to accommodate simultaneously the new 
CSB result of our study and that of a recent cross-section 
measurement of the isospin forbidden reaction dd — > air 
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